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ABSTRACT 



Context. Lithium abundances in open clusters are a very effective probe of mixing processes, and their study can help us to understand 
the large depletion of lithium that occurs in the Sun. Owing to its age and metallicity, the open cluster M67 is especially interesting 
on this respect. Many studies of lithium abundances in M67 have been performed, but a homogeneous global analysis of lithium in 
stars from subsolar masses and extending to the most massive members, has yet to be accomplished for a large sample based on 
high-quality spectra. 

Aims. We test our non-standard models, which were calibrated using the Sun with observational data. 

Methods. We collect literature data to analyze, for the first time in a homogeneous way, the non-local thermal equilibrium lithium 
abundances of all observed single stars in M67 more massive than ~ 0.9 Mq. Our grid of evolutionary models is computed assuming 
a non-standard mixing at metallicity [Fe/H] = 0.01, using the Toulouse-Geneva evolution code. Our analysis starts from the entrance 
in the zero-age main- sequence. 

Results. Lithium in M67 is a tight function of mass for stars more massive than the Sun, apart from a few outliers. A plateau in 
lithium abundances is observed for turn-off stars. Both less massive {M < 1.10 Mq) and more massive (M > 1.28 Mq) stars are more 
depleted than those in the plateau. There is a significant scatter in lithium abundances for any given mass M < 1.1 Mq. 
Conclusions. Our models qualitatively reproduce most of the features described above, although the predicted depletion of lithium is 
0.45 dex smaller than observed for masses in the plateau region, i.e. between 1.1 and 1.28 solar masses. More work is clearly needed 
to accurately reproduce the observations. Despite hints that chromospheric activity and rotation play a role in lithium depletion, no 
firm conclusion can be drawn with the presently available data. 

Key words. Stars: fundamental parameters - Stars: abundances - Stars: evolution - Stars: interiors - Stars: solar-type 



1. Introduction 

Lithium is destroyed by proton capture at temperatures above 
2.4 X 10^ K. In Sun-like stars, material in layers at this temper- 
ature can be mixed with photospheric layers through convection 
and diverse processes such as diffusion, me ridional circulation, 
and internal gravity waves (e.g. lTalonll2008l) . As a consequence, 
stars undergoing these processes have lower lithium abundances 
than their initial values. This phenomenon is known as lithium 
depletion, and makes lithium abundance one of the most effec- 
tive probes of mixing processes in stars. 

Astronomers have long dedicated huge observational and 
theoretical efforts to understand how and when lithium is 
depleted in stars. Depletion of lithium can occur during both 
the pre-main sequence (PMS) and the main sequence (MS 
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lifetim e of solar-type stars (e.g. D'Antona& Mazzitelli"198 
1994tlD^ Antona & Montalban 2003; Boesgaard & Budge 1988£ 



iMontalban & SchatzmanI 119961: IChen et all 12001 



Piau & Turck-Chi ezel l2002t 



Lambert & Reddy, 2004^ 



"Sestito & RandichI 120051: ISestito et alj l2006l) 



"'Theado & Vauclair ^2003b 
Takeda & Kawanomoto 2005 ; 
Lithium de- 
pletion by rotationally driven mixing has been related to the 
rotational history of stars, which, in turn, has been related to 
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planet formation (e.g. lBouvierll2008l) . The influence of planets 
is the subject of much debate. There are claims of stronger 
lithium depletion in planet-hos ting stars but only in a very 
narro w range around solar T^fi dGonzalezI 120081 : llsraelian et al.l 
l2009h . However, This evidence has been challenged by other 
authors, who claim that no difference in lithium abundance 
is detected between stars with and wi thout detected planets 
when unbiased samples are analyzed ([Melendez et all 120101: 
iGhezzi et al.l 120101 : iBaumann et alJ l2010l) . Furthermore, at a 
given mass and metallicity there is a clear relation between 
lithium abundances and age, independently of the star being a 
planet host or not (Melendez et al. 2010: Baumann et al. 2010). 

Obviously, there are two main parameters favoring lithium 
depletion: the depth, hence the temperature to which the mixing 
extends and the length of time for which the process is active, 
i.e. the age of the star. However, the dependence of the former 
parameter on stellar mass, age, and metallicity is not simple, nor 
is there a general consensus on which additional parameters are 
involved. Despite decades of observations in clusters and field 
stars, no clear and widely accepted picture has yet been drawn. 

Lithium abundance in open clusters is a precious tool for 
understanding mixing processes because stars in open clusters 
share the same age and initial chemical composition, therefore 
allow us to see how lithium depletion varies with stellar mass and 
to test non-standard models of stellar mixing. A crucial question, 
on which observations of clusters can cast light, is whether there 
is a spread in lithium abundance for stars with the same fun- 
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damental parameters, i.e. mass, chemical composition, and age, 
which is to say, whether lithium abundance in a given open clus- 
ter is a tight function of the stellar mass. In case it is not, careful 
study of the deviating stars would be important to constrain the 
parameters that cause the spread of lithiu m abundances, such 
as for example the initial rotation velocity (ICharbonnel & TalonI 
I2005h . 

Some early observations indicated that there is a certain 
amount of spread in the lithium abundances versus mass dis- 
tribution in a few clusters, and its absence, or at least the lack 
of evidence for it, in many others. Example s of clusters be- 
longi n g to the former group a re the Pleiades (Duncan & Jones 
1983t fSoderblom et al.lll993ah and Praesepe (Soderblomet al 



ICastroetall (1201 
(l2007h . iJones et al 



) iPasQuini et al.1 (l2008h. [Randich et aH 
(I1999I) . and iBalachandraiil (Il995h . When 
multiple measurements were available for a single object, 
we adopted the most recent. Some of the data i nclud ed in 
our compil ation, namel y data from ICastro et al.l (1201 lb and 
Table 2 in iJones et all (Il999l) . are refinements of previous 
abundance determinations. The final sample amounts to 103 
stars. For the solar twin YBP 1194, Onehag et al.l (1201 lb used 
the highest quality spectra available to date to determine its 
lithium abu ndance, although th eir result is essentially identical 
with that of lCastro et al.l (1201 lb . which is based on lower quality 
spectra; thus, either of the sources can be adopt ed to measure 
the lithium abundance of YBP 1 194. The data in lPasauini et alJ 



1993bl) for low-mass stars (M < 0.9 Mq). However, many of (Il997h. iGarcia Lopez et al.l ( Il988h . ISpite et al 



the early claims about a spread in the lithium ver sus mass distri 
bution in ope n clusters have lately been revised. IXiong & Deng! 
(I2005L l200(5h found that most, if not all, of the lithium disper- 
sion for a given mass in members of o-tmb Per and Pleiades, 
can be accounted for by inhomogeneous reddening, stellar spots, 
and stell ar surface activity, rather than by a genuine variation. 
iKing et a l. (2000), analyzing Praesepe stars, ascribed most of 
the scatter to variations in activity-regulated ionization of the 
lithium atom. In support of this hypothesis, they mentioned the 
lesser amount of scatter in t wo older and chromos pherically 
quieter open clusters: Hyades (iThorburn et al.l[l993l) and M34 
dJones et al.lll997i) . Now, we can add to the list several cluster s 
older tha n the Hyades, such a s NGC752 dSestito et al.l |2004. 
NGC 188 (iRandich et alJl2003h . and Berkelev32 (iRandich et al] 
|2009|) . There are, however, examples of old clusters presenting a 
significant amount of scatter One of them is NGC3680, but the 
scatter i n this cluster is especially at t emperatures higher than ~ 
6500 K ( lAnthonv-Twarog et al.l2009h . and too few stars at lower 
temperature have been analyzed. Another old cluster apparently 
showing a spread in the lithium abundance distribution is M67, 
which is a well-studied cluster, subject to a number of lithium 
abundance investigations. However, a comprehensive study in 
which lithium abundances were analyzed for a large number of 
members in all evolutionary stages, taking advantage of already 
available high signal-to-noise ratio (S/N) spectra, is still lacking. 
Since M67 has an age of about 3.9 Gyr (e.g. 
IVandenBerg et al.l 120071: ICastro et al.l 1201 II) a nd about so 



lar metallicity (e . g. [ Tautvaisiene et al. l2000t IRandich et al 



lar metallicity (e . g. [ iautvaisiene et aL UO{J( 
20061: iPace et al.l 120081: iPasquini et al.l l200i 



Onehag et al. 



20111) . it can provide insight into the behavior of lithium in 



relatively old solar type stars, which is important to assess 
why some old solar a nalogs in the field see m to have too much 
lithium for their ages (iBaumann et al.ll2010l) . 

In this work, we homogenize all previous studies on M67 
stars by obtaining a consistent set of temperatures, applying cor- 
rections to the lithium abundances when needed (due to revised 
effective temperatures), and taking into account non-local ther- 
mal equilibrium (non-LTE) effects on the lithium abundances. 
We analyze stars from the MS to the red giant branch (RGB), 
comparing our models to the observed behavior of lithium abun- 
dance as a function of mass. In Sect. |2] we describe the data. 
In Sect.|3j we compare the effective temperatures from different 
references used in the paper In Sect.|4] we present details of our 
stellar models, and in Sect. |5]we discuss our results. Finally, in 
Sect. |6] we give our conclusions. 

2. Sample and data 

Our database is a compilation of Uterature sources of lithium 
abundance measurements, namely ICanto Martins et aP (1201 ll) . 



r-i— I (II987h . and 

iHobbs & Pilachowskil ( Il986l) were initially considered for 
the present compilation but eventually not used, since their 
whole sample was later reanalyzed in some of the works cited 
above. We also did no t include in our compilation data from 
iDelivannis et al.l (1 19941) . which is a study of tidally locked bi- 
naries, whose mechanism of lithium depletion may completely 
differ from that of single stars. 

In some works, no quantitative estimations of the error in the 
lithiu m abundances were gi ven. All these cases are discussed 
below. iPasquini et alJ (l2008l) stated that all the stars with upper 
limits in their sample may have a lithium abundance compa- 
rable to the solar one. Following this piece of information, we 
raised their upper limits (which in some cases were as low as 0.5 
dex) to 1.0 dex, which is more reasonable considerin g the qual- 
ity of the spectra available to IPasquini et al.l (l2008l) . However, 
this correction hardly ch anges the gene r al pic ture drawn us- 
ing data in Table A. 2 of IPasquini et al.l (|2008|) . which points 
to a significant spread in lithium abundances at about one so- 
lar mass. Their estimation of the uncertainty in lithium abun- 
dance due to the uncertainty in the equivalent width, for a star 
with A(Li) = 2.2, is ±0.04 dex, and no evaluation is given for 
stars with lower abundances. Since lithium abundance is ex- 
tremely sensitive to temperature, the uncertainty in the latter may 
contribute to that of the former in a non-negligible way, even 
when it is extremely low. This contributio n was not quantified by 
IPasquini et al. (2008). Castro et al.l( l201 II) . in reanalyzing the so- 
lar twins of .Pasquini et al.. (20081) with spectral synthesis, found 
a global error of 0.1 dex. Their stars have Li abun dances that are 
lower than the average of the sample studied bv IPasquini et al.l 
(|2008|) . We conservatively a ssumed this error (0.1 dex) for all 
sta rs inlPasQuini et al.l (l2008l) . 

I Jones et al] d 1 9991) estimated an error of 0.05 dex in lithium 
abundance due to the uncertainty in the equivalent width, and 
an equal error due to the uncertainty in the temperature of 50 K. 
We were conservative once again, and we summed linearly these 
two sources of error and assumed for all lithium abundances an 
uncertainty of 0. 1 dex. This same value was also assumed for 
their measurements made by reanalyzing literature data. 

Information on photometry and membership was taken, 
w henever possible, i .e. for 95 out of 103 stars, from the work 
of lYadav et al.l (l2008h . otherwise it was taken from the reference 
of the stud y on lith ium abundance. Of the 95 stars of our sample 
studied bv lYadav e t al. (2008), 87 had radial velocity from the 
same source. For the remaining 8 stars, Yadav et al.'s member- 
ship information was based on photometry and pro per motions . 
Six of the 8 stars of the final sample not studied bv lYadavet al.l 
(l2008l) . have radial velocity data available from the lithium abun- 
dance study, therefore only 10 stars do not have radial velocities. 
We discarded stars with radial velocities outside the range of val- 
ues from 30 to 38 km s"' and removed from our sample the blue 
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Fig. 1. Color-magnitude diagram for M67 superposed upon the 
isochrone, which corresponds to an age of 3.87 Gyr The red 
filled circles represent our work sample, i.e. stars with lithium 
abundance measure ments. Stars with photometry available from 
lYadav et all (l2008h . but with no lithium abundances, are dis- 
played, for comparison with the isochrone, as black points. The 
squares represent the deviant stars. The roman numbers I, II, and 
III are described in Sect. 5.1 and correspond to the temperatures 
of, respectively, 6095, 6123, and 4966 K. 



stragglers Sanders 2204 and Sanders 997 and other known bina- 
ries. 

In Fig. [1] we plot a color-magnitude diagram of M67 stars, 
and the isochrone calc ulated from our mod els, with the param- 
eters for M67 found in lCastro et aTl(l20Ill) . i.e. [Fe/H] = 0.01, 
(m-M) = 9.68, and E{B-V) = 0.02, ar idageof 3.87 Gyr. B lack 
dots are the stars from the sample of lYadav et aP (l2008i) . red 
filled circles are the stars studied in this paper, and open squares 
are deviant stars (see Sect. 15. 3) . Roman numbers I, II, and III are 
defined in Sect. |5] 

Photometric and membership information, are available at 
CDS in electronic form (Table 1). 



3. Temperatures and corrected lithium abundances 

Lithium abundances are strongly sensitive to the effective tem- 
perature. To check for the dependence on the adopted effective 
temperatures, we compared different estimates, using the litera- 
ture source from which the lithium abundances were also taken, 
a photometric calibration, and the isochrone. To determine the 
effective temperature by means of the isochrone, we proceeded 
in the following way: we plotted in the same CMD data points 
for each star and the isochrone computed for M67 age and metal- 
licity, shifted according to its distance modulus and reddening. 
For each stellar data point, we then took the closest point on the 
isochrone, and adopted the corresponding physical parameters, 
i.e., in the case under consideration, effective temperature. For 
the remainder of the paper, we refer to this kind of temperature 
estimation as either "isochrone temperature", or Tiso with the 
further specification of the isochrone used. "Isochrone masses" 



were computed in the very same way. Summarizing, we have the 
following four temperature estimations: 

1 . the Teff adopted in the paper from which the lithium abun- 
dance was taken (Torig); 

2. the isochron e temperature adop ting models and parameters 
for M67 as inlYadav et al.l(l2008l) . who in turn adopted BaSTI 
models dPietrinfemi et al.ll2004h (risoBaSTi); 

3. the isochrone temperature adopting isochrones and parame- 
ters for M67 in Castro et al. (20 11) (riso-castrol; 

4. using a photometric calibration, from ICasagran de et al 



201C ) for dwarfs and sub-giants, and from IKucinskas et al 



( 120051) for giants {Tcniib or "calibration temperature"). 



A comparison between the different temperature evaluations 
shows that the differences are significant. When comparing Torfg 
with the other temperature evaluations, the heterogeneity of 
the data, especially in the assumed reddening for M67, is cer- 
tainly a main cause of these differences, and one of the aims 
of the present study is to homogenize the different determina- 
tions. The mean differences between Tgrig on the one hand, and 
T'iso-BaSTi, T'iso-Castro, and Tcaiib on the other, are, respectively, 19, 
50, and 95 K. The differences between calibration and isochrone 
temperatures arise, instead, from the isochrone in the CMD leav- 
ing most of the points on the right (cooler) side. This is because 
the calibration temperature is affected by the presence of un- 
detected companions, which move the data points brighter and 
cooler in the CMD. In computing isochrone temperatures, we 
exploit the fact that the stars belong to the same cluster The 
isochrone on the CMD is the locus where we would expect all 
the photometric data points to be if there were no photometric 
errors and no multiple systems. Ti^o is the temperature corre- 
sponding to this locus, therefore it should be unaffected by the 
presence of undetected companions. This is confirmed by the 
overall go od agreement betwe e n Tj^g and Tone for t he 82 stars 
studied inlCanto Martins et all (l201lh . ICastro et al.l (l201ll) . and 
iPasquini et al.l (l2008b . In these stars, T^rig is based on spectro- 



scopic analysis and is therefore very precise. We label Tspec this 
subset of Torig based on spectroscopic analysis. The linear best 
fit to the data points in the Tspec versus isochrone temperature 
graph is very close to the identity: Tiso-BaSTi - 0.99 ■ Tspec - 2.5 
K. We can confidently claim that the errors in the isochrone tem- 
peratures are random in nature. 

We note that the two different evaluations of effective tem- 
perature using either isochrone, i.e. Tiso-BaSTi and Tiso-castro, are 
consistent with each other within the margins of error for most 
of the stars. The only significant differences arise for sub-giant 
stars and, to a lesser extent, for turn-off stars. 

Corrections to lithium abundance measurements were ap- 
plied according to the differences between Torfg and, respec- 
tively, Tiso-BaSTi, Tiso-castro, and Tcaiib- At the same time, we also 
took into account NLTE effects usin g the grid of NLTE correc- 
tions computed bv lLind et alJ (l2009h . 

All temperature estimations, corrected and uncorrected 
lithium abundances, along with the photometric and member- 
ship information, are available on-line in Table 1 . 

In Fig. 121 we plotted the lithium abundances as a function of 
mass for both the stellar data points and two models. Red filled 
circles with error bars are the stars of our sample with detected 
lithium line, and green filled triangles represent upper limits. 
Open squares with error bars and, only in the case of one up- 
per limit, the open triangle, represent the deviant stars discussed 
in Sect. 15.31 The models are discussed extensively in Sec t.lH 

As me n tioned above^ [Canto Martins et al] ( 1201 lb . 

ICastro et aP (1201 ih . and iPasquini et al., (.2008;) computed 
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Fig. 2. Lithium abundances as a function of stellar mass deter- 
mined from the isochrone. Red filled circles represent stars with 
lithium abundance redetermined. The green triangles represent 
the lithium abundance upper limits. The open squares represents 
the deviant stars as discussed in Sect. 5.3. The continuous line 
corresponds to the lithium abundance predicted by our models 
with rotation-induced mixing calibrated for the Sun at the age of 
M67. The dashed line represents the same model shifted by 0.45 
dex to match the observed lithium plateau. The dotted line rep- 
resents models with rotation-induced mixing calibrated on the 
Hyades. The roman numbers I, II, and III are described in Sect. 
5.L 



temperature by carrying out a detailed spectroscopic analysis. 
Their temperature evaluation is more precise than that obtained 
using photometric information. Therefore, for the 82 stars 
studied in these works, we used the values given in the literature 
to produce Fig. |2] In contrast, for the remaining 21 stars, we 
used riso-Castio, and the corresponding corrected value of the 
lithium abundance. For all of the stars, the mass adopted is the 
isochrone mass, which was comp uted using the isoch rone and 
the parameters forM67 adopted in lCastro et al.l(l201lb . 



4. Stellar evolutionary models 



For this study, stellar evolutionary 
using the Toulouse-Geneva stellar 



(iHui-Bon-Hoal 



2008) 



models were computed 

evolution code TGEC 

Details on the physics of these models 



can be found in Richard et al. ("19 961120041) . Ido Nascimento et al.l 
(120001) . and lHui-Bon-Hoa ( 2008i) 



4.1. Input physics 

We used the OPAL2 001 equation of state by 
Rogers & Navfonov (l2002h and the radiative opacities by 
Iglesias & RogersI (Il996l) . completed wi th the low temperature 
atomic and molecular opacities by [Alexander & FergusonI 
(Il994l) . The nuclear reaction s were described by the analytical 
formulae of the NACRE (lAngulo et al.l 1 19991) compilation, 
taking into account the three pp chains and the CNO tricycle 



with the iBahcall & Pinsonneauld (Il992[) screening routine . 
Convection was treated according to the iBohm-Vitensel (Il958h 
formalism of the mixing length theory with a mixing length 
parameter a - l/Hp - 1 .74, where I is the mixing length and 
//p the pressure height scale. For the atmosphere, we used a 
gray atmosphere following the Eddington re lation, which is a 
good approximation for MS solar-type stars (IVandenBerg et al.l 
2008 ). 

The abundance variations in the following chemical species 
were computed individually in the stellar evolution code: H, He, 
C, N, O, Ne, and Mg. Both Li and Be were treated separately 
only as a fraction of the initial abundance. The heavier elements 
were gat hered in a mean species Z . The initial composition fol- 
lows the iGrevesse & Noelsl (Il993h mixture with an initial he- 
lium abundance Fini = 0.268. We chose to use the "old" abun- 
dancesofjGrevesse & Noelsl (1 1993h instead of the "new" mixture 
of lAsplund et al.l ( 120091) . This choice was motivated by the dis- 
agreement between models computed with the new abundances 
and the helioseismic inversions for the sound-speed profile, the 
surface helium abundance, and the convective zone de pth (e.g. 
Serene lli et al. 2009). Furthermore, according to Caff au et aP 
(l2009h . the solar metallicity using their own three-dimensional 
analysis is given by the value s Z - 0.0156 and Z/X = 0.0213, 
which are closer to those of IGrevesse & Noelsl (Il993l) . In all 
cases, the solar abundances remain uncer tain. We note that 
the accretion of metal- poor material (e.g. ICastro et al.l 120071 : 
iGuzik & Mussacldl2010l) as suppor ted by the lack of refr actory 
elements in the solar atmosphere (iMelendez et al.l l2009l) . may 
help to reduce the discrepancy between the solar model and he- 
lioseismic data. Thus, the low solar abundances of lAsplund et al.l 
(i2009i) may actually be representative of the solar photosphere, 
while in the solar interior the abundances may be higher 
Diffusion and rotation-induced mixing. 

The microscopic diffusion is computed with the atom test 
approximation. All models include gra vitational settli n g wit h 
diffusion coefficients computed as in iPaquette et al.l (Il986l) . 
Radiative accelerations are not computed here, since we fo- 
cus only on solar-type stars where th eir effects are small 
when mixing is taken into account dTurcotte et al.l 119981: 



iDelahave & Pinsonneaultl 20051). Rotation-induced mixing is 



computed as described in lTheado & Vauclaij ( 2003al). This pre- 
scription is an e xtension of the approach of IZahnI (1 19921) and 
iMaeder & ZahnI (1199 81 and introduces the feedback effect of the 
ju-currents in the meridional circulation, caused by the diffusion- 
induced molecular weight gradients. It introduces two free pa- 
rameters in the computat ions: Ch and aturb (cf. Eq. (20) of 
iTheado & Vauclain 12003 al) . The ev olution in the rota tion pro- 
file follows the Skumanich's law (ISkumanichI 1 1 972h with an 
initial surface rotation velocity on the zero-age main-sequence 
(ZAMS) equal to Vi = 100 km.s ', which roughly corre- 
sponds to the mean rotation velocity of stars hotter than about 
7000 K in t he Hyades according to the statistical study of 
iGaigd (Il993l) . The parameters of the Skumanich's law are cal- 
ibrated to match the solar rotation velocity at the solar age 
(~ 2 km.s '). Other prescriptions were tested by other authors 
to model the lithium de struction. |Charbonnel & Talon (19991) 
and lPalacios et alJ (l2003l) included angular momentum transport 
induced by mixing. However, since a rotation-induced mixing 
alone cannot account for the flat rotation profile inside the Sun, 
these authors later introduced the possible effect of internal grav- 
ity waves triggered at the b ottom of the convective zone (see e.g. 
iTalon & Charbonnelll2005l) . which allows the hot side of the Li- 
dip to be reproduced. Other authors suggested that the internal 



Pace et al.: Lithium in M67: from the MS to the RGB 



magnetic field is more impo rtant than internal waves in trans- 
porting angular momentum (iGough & Mclntvrall998h . In any 
case, when applied to the solar case, all of these prescriptions are 
able to reproduce the lithium depletion obs erved in the Hyades, 
and the results are ultimately quite similar (iTalon & Charbonnell 
[T99ilTtiiado & Vauclaidl2003bi) . 

We also include a shear layer belo w the c onvective zone, 
which is treated as a tachocline (see ISpiegel & Zahn 1992) . 
This layer is parametrized with an effect ive difiiision coefficient 
that decreases expon entially downwards (iBrun et al.ll998[ri999t 
iRichard et aLll2004l) : 



Aacho = £>bczexp(ln2 — -^j 



where Dbcz and rbcz are the value of Dtacho at the bottom of the 
convective zone and the radius at this location, respectively, and 
A is the half width of the tachocline. Both Dbcz and A are free pa- 
rameters and the absolute size of the tachocline (i.e., A/7?» where 
Rt is the radius of the star) is supposed to be constant during the 
evolution. An overshooting of parameter ofov - 0.01//p has been 
included in the models that develop a convective core. 

4.2. Models and calibration 

We calibrated a model of 1 .00 Mq to match the observed solar ef- 
fective temperature and luminosity at the s olar age. The calibra- 
tion method of the models is based on the iRichard et al.l ( 1 19961) 
prescription: for a 1.00 Mq star, we adjusted the mixing-length 
parameter a and the initial helium abundance ¥[„■, to reproduce 
the observed solar luminosity and radi us at the solar age. Th e ob- 
served values that we used are those of lRichard et al.l(l2004l) . i.e., 
Lq = 3.8515 + 0.0055 x lO^"* erg.s"', R^ = 6.95749 + 0.00241 x 
10'° cm, and t^ - 4.57 + 0.02 Gyr For the best-fit solar model, 
we obtained L = 3.8501 x 10^^ erg.s"' and R = 6.95524 x 10'° 
cm at an age f = 4.576 Gyr. 

The free parameters of the rotation-induced mixing deter- 
mine the efficiency of the turbulent motions. They are ad- 
justed to produce a mixing that is both: 1) efficient and deep 
enough to smooth the diffusion-induced helium gradient below 
the surface convective zone, thus improving the agreement be- 
tween the model and seismic sound speed profiles; 2) weak 
and shallow enoug h to av oid the destruction of Be. Following 
iGrevesse & Sauvall (Il998h . the Be abundance of the Sun is 
A(Be) - 1.40 + 0.09. We obtained a slight Be destruction by 
a factor of 1.33 with respect to the meteoritic value, which is 
well within the error in the solar Be abundance. 

The calibration of the tachocline allows us to reach the solar 
hthiu m depletion (A(Li) - 1.10 + 0.10 (e.g. IGrevesse & Sauvall 
Il998h and for our best-fit solar model we obtained A(Li) - 1 .04. 
We also checked that the sound velocity profile of our best-fit 
model is co nsistent with that d educed from helioseismology in- 
versions by iBasu et aD (1 19971) . Our calibration is in excellent 
agreement with helioseismology, more accurately than 1% for 
most of the star, except in the deep interior, where the discrep- 
ancy reaches 1.5%. 

We computed a grid of evolutionary models of masses in the 
0.90 to 1.34 Mo range, with a step of 0.01 M© in mass, and a 
metallicity of [Fe/H] = 0.01, which is a si mple average of differ - 
ent estimates of the metallicity of M67 (see lPasquini et al.l2008l) . 
We ran the models from the ZAMS to the top of the RGB for the 
most massive stars. The input parameters for all the models are 
the same as for the 1.00 M© model. 

There is no a priori reason why the calibration of the pa- 
rameters of the extra mixing for the model of 1 .00 M© should 



also hold for different masses, and at di fferent times during 
the evo lution of the stars. In their work, iTheado & Vauclain 
(l2003b|) adjusted the free parameters of the rotation-induced 
mixing with the TGEC for each one of the models of differ- 
ent masses to obtain the correct lithium depletion at the age 
of the Hyades. They noted that the horizontal diffusion coeffi- 
cient is very mass-dependent. The uncertainty associated with 
these parameters led us to analyze their impact on the de- 
struction of lithium as a function of stellar mass. To do so, a 
first set of models was calculated to reproduce the profile of 
lithium abundance as a function of mass for the Hyades open 
cluster. Our Hyades sample is a compilation of EW measure- 
ments of the lithium fine at 6708 AA from several sources 
(iRandich et al J 120071: lThorburnetal.1 119931: ISoderblom et all 



19901: iBoesgaard & Budgd 119881: iBoesgai^d & Tripiccol 1198 



i 



Duncan & Joneslll983l:lRebolo & Beckmarjll988t) . All the stars 
identified as binaries in iThorburn et al.l (1 19931) were excluded. 
When more than one reference was available for a star, we used 
the most re cent. Masses and temp eratures were taken either di- 
rectly from iBalachandranI (1 19951) , or computed from the rela- 
tionship between B-V and, respectively, mass and temperature, 
obtained from the data in the same paper The result of the com- 
pilation is available at CDS in electronic form (Table 2). The 
parameters Ch and o-tuib of the rotation-induced mixing of the 
models, with masses from 0.800 to 1.400 M©, with a step of 
0.050 Mo, are calibrated to reproduce at the age of the Hyades 
(625 Myr) the observed destruction of lithium. In these calibra- 
tions, we did not include the tachocline that is used in the code to 
calibrate solar models, because in the solar mass stars the bottom 
of the convective zone is very close to the lithium destruction 
layer (JTheado & Vauclair 2003b) . The same calibration is then 
used to calculate a set of models of M67 stars. The values of the 
parameters Ch and o-tuib for all calibrations are given in Table |3] 

5. Discussion 

The observed spread in the lithium abundance versus mass distri- 
bution of M67 for stars of about one solar mass or lower, which 
can be clearly seen in Fig. |2] suggests that an extra variable (in 
addition to mass, age, and metallicity) is at work, leading to dif- 
ferent lithium depletion in cluster stars of the same mass. The 
fundamental question here (since the spread is present) is which 
other variable may cause a large range of lithium abundance at a 
fixed mass? 

Our database includes stars ranging from the MS to the sub- 
giant and giant branches, so is well-suited to both address the 
issue of spread and test models of non-standard mixing for dif- 
ferent masses. 

Models that point to rot ational mixing as th e major cause 
of lithium depletion (see e.g. |Pinsonneaultll20 1 Q. for a review), 
predict that this only depends on the age, metallicity, mass, and 
the initial angular momentum of the star This last parameter, 
however, may vary for otherwise similar stars, introducing a 
spread in the lithium abundances. The variation in the initial ro- 
tation velocity by a factor two in our models has a weak effect 
of about 0.05 dex of magnitude for the lithium destruction, be- 
cause the Skumanich's law used implies that there is a strong 
drop in the rotation velocity very early on in stellar evolution. 
The Skumanich's law is empirical, and we should in the future 
include in the code the transport of angular momentum to ensure 
a more reliable estimation of the spread of the lithium abundance 
due to a possible spread in the initial angular momentum. 

The depth of the surface convection zone and the nu- 
clear burning depend strongly on the total stellar mass 
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Table 3. Parameters of the rotation-induced mixing in TGEC models. 



Calibration 


Mass (Mo) 


Ch 


Q'turb 


Sun 


1.000 


9000 


1.00 




0.800 


3400 


4.00 




0.850 


10000 


5.00 




0.900 


10000 


5.50 




0.950 


10000 


5.60 




1.000 


10000 


5.15 




1.050 


7000 


4.25 


Hyades 


1.100 


5000 


3.30 




1.150 


4500 


2.50 




1.200 


2500 


2.10 




1.250 


2546 


2.15 




1.300 


2600 


2.75 




1.350 


6500 


3.50 




1.400 


10000 


4.50 



ddo Nascimento et al]|2009l) . thus the most appropriate way to 
investigate the abundance of lithium in these M67 stars is first 
to divide our sample of stars in three different groups of mass 
range. MS stars with M < 1.1 M©, stars close to the turnoff with 
1.1 Mo< M < 1.28 Mo, and evolved sub-giant and giant stars 
withM>1.28Mo. 



5. 1 . Lithium abundance and mass 

In Fig. |2] we compare our model predictions with our inferred 
stellar masses and lithium abundances for our sample of M67 
stars. The solid line represents an isochrone constructed with 
all the models including the effects of atomic diffusion, gravi- 
tational settling, and rotation-induced mixing calibrated for the 
Sun as described in Sect. [H T he ag e of the models is 3.87 Gyr 
as determined infCastro et aLJ (1201 lb . Initial lithium abundances 
were chosen to equal A(Li) = 3.26, the estimat ed initial solar 
value based on meteorites dAsplund et al.l 120091) . This lithium 
abundance is 0.05 dex lower than the previo usly accepted mete- 
oritic lithium abundance (A(Li) = 3.31, see lGrevesse & Sauvall 
Il998l) . The dashed line represents an isochrone for the same 
models but where the initial lithium abundance was rescaled by 
-0.45 dex to fit the plateau at A(Li) = 2.5. This suggests that 
the initial M67 lithium abundance may have been lower than 
3.26, or that our models are not depletin g enoug h lithium. We 
note, however, that according to Delivan nis et al.l ([l994), based 
on a short-period tidally locked binary in M67, the initial lithium 
abundance in M67 was at least A(Li) = 3.0. This would reduce 
the shift necessary to match our models with the lithium content 
in the plateau to only about -0.2 dex. However, the hypothesis 
of a lower initial lithium content seems weak and unrealistic. We 
are left with an insufficient destruction of lithium by rot ationally 
driven mixing mechanisms. iTheado & Vauclaid (l2003bl) showed 
that for masses higher than 1 .0 Mq, it is necessary to increase 
the values of these parameters as a function of mass to account 
for the lithium depletion at the age of the Hyades. To analyze 
the influence of these parameters on the lithium destruction, we 
calculated another set of models, in which the rotation-induced 
mixing parameters was calibrated for each mass, as described in 
Sect.|4] The isochrone at the age of M67 of these models is repre- 
sented by the dotted line in Fig.|2l These models reproduce more 
closely in a quantitative way the profile of lithium destruction in 
M67, but the arbitrary calibration of the mixing parameters for 
each mass is a very unsatisfactory method. There appears to be 
no clear relation between the mixing parameters and mass that 
might have a physical significance. 



For stars with masses lower than 1.1 Mq (roman num- 
ber I in Fig. [T] and Fig. |2]i, the lithium depletion progres- 
sively increases toward lower masses. A dispersion is ob- 
served in this mass r ange and confirm e d by several authors 
(Pasquini et all 1 19971: lJonesetal.1 Il999t iRandich et all l2002l 



I2007t IPasquini et al 
upper limits used 



2008h . We note that in this region, some 
may be con sidered as too optimistic. In 



this context, lOnehag et al.l (1201 ih analyzed the M67 solar twin 
YBP 1194 using spectra of good quahty (R ~ 50, 000, S/N = 
160), but even with such a spectrum they found that it is chal- 
lenging to robustly determine a lithium abundance in the solar 
twin YBP 1 194. A careful analysis of solar analogs in M67 us- 
ing high quality spectra is urgently needed for stars around 1 
Mq. However, we note that four stars at 1.01 M© span as wide 
a range of lithium abundances as 1 . 1 dex, and one of them is 
a revised upper limit. A considerable amount of scatter in the 
data in this region remains after raising the upper limits, whose 
real presence should be considered highly probable based on the 
presently available data. Stellar rotation and the history of angu- 
lar momentum induce it either directly through rotational mixing 
or indirectly by drivi ng other processes such as diffus ion and in- 
ternal gravity waves (iMontalban & Schatzmanll2000l) . Planetary 
accretion may also be at work in sol ar-type stars, and i nduce 
mass-independent lithium destruction dTheado et al.ll2010l) . 

The low lithium abundance observed in low-mass stars in 
M67 also exists because these stars have had more time to burn 
their lithium during the PMS and have a surface convection zone 
that both retreats more slowly and ends up at greater depths on 
the MS. We note that, for this range of masses and these evo- 
lutionary stages, the standard models predict that the bottom of 
the convective zone is not hot enough to account for the decrease 
in lithium, even though the observed low lithium abundances 
clearly indicates the need for a more realistic representation of 
transport mechanisms in low mass stars. 

Stars with masses 1.1 Mq< M < 1.28 Mq (between ro- 
man number I and II in Figs. [T] and |2| are close to the turn off 
(or about to leave the MS) and at these masses the surface con- 
vection zone and the mixed layers below the convection zone are 
too thin to support temperatures high enough to cause the nuclear 
destruction of lithium. Furthermore in this mass range the PMS 
lithium depletion tapers off to essentially zero. We can therefore 
easily see why there should be a plateau with a nearly constant 
lithium abundance for this mass range. One important test of our 
model will be in fact to reproduce the general morphology of the 
lithium depletion behavior 
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This is the region where the offset between our models and 
observations becomes apparent, as discussed above. 

Stars with masses above 1.28 Mq are sub-giants follow- 
ing an evolutionary path from the turn-off to the RGB. In Figs. 
[T] and 12] this stage is marked by the roman numbers II and III. 
The lithium abundance predicted for standard mod els of sub- 
giant and giant sta rs is rnainly controlled by dilution (llbenll 19661 
Il967a.b: Scalo & Milledl980l) . This process was described orig- 
inally by llbeni (Il965h and appears when a star evolves off the 
MS. When the stars crosses the sub-giant branch and climbs 
the RGB, the surface convective zone increases its fraction in 
mass. Lithium-poor material rises to the surface and is mixed 
with Li-rich material. This process stops when the convective 
zone reaches its maximum size. From our models, the pre- 
dicted value for lithium-abundance post-dilution is A(Li) = LO 
dex. The observations clearly show the need for a non-standard 
transport processes for su b-giants and evolved stars of M67. 
ICanto Martins et all ( 1201 lb discussed the depletion of the sur- 
face Uthium abundance and whether there is an additional non- 
standard transport process related to the transport of matter and 
angular mo mentum by meridional cir culation and shear-induced 
turbulence (iRvan & Delivannislll995h . 



5.2. Prediction by other models 

Our models successfully reproduce qualitatively the lithium de- 
pletion seen in different mass regimes in M67 stars, as shown 
in Fig. 121 In Fig. [3] the evolution with age at different masses is 
shown. We note that f ully independent non-standard models by 
IXiong & Dend (l2009l) . also show similar trends between lithium 
depletion, mass, and age. 

As shown in Fig. 1 of IXiong & Dend ( |2009|) . at 4 Gyr, there 
is a large depletion in the lithium for models of 1 Mq. As in our 
models, the less depleted stars are those with masses som ewhat 
higher than LI Mq, although the IXiong & Dend (l2009l) model 
also fails to predict exactly the peak or plateau in lithium abun- 
dance for M67 stars close to the turnoff, requiring probably an 
additional lithium depletion of 0.3 dex. Stars with masses higher 
than 1.3 M© show severe depletion at 4 Gyr, probably because 
they should be in th e sub-giant or giant ph ase. Thus, qualita- 
tively, the models by IXiong & Dend ( l2009l) also reproduce the 
M67 observations. 



5.3. Deviant stars 

Although overall lithium abundance appears to be a tight func- 
tion of mass for stars more massive than 1 .06 M©, a few objects 
dramatically depart from the general trend. 

These objects are marked by open square symbols in Figs. [1] 
and|2] and are discussed below. 

YBP 1075, YBP 750, YBP 769: These objects, despite be- 
ing quite distant from the isochrone in the CMD diagram (more 
than 0. 1 mag, limit under which 80% of the sample lies), are 
all secure members (99 or 100% of probability, based also on 
radial velocity). They are therefore likely to have companions 
too small to be detected, but bright enough to affect color and 
brightness. Their lithium abundance deficiency, with respect to 
the mean trend with the mass, is about 2.3, 1.0, and 0.5, dex 
respectively, much more than errors can possibly explain, and 
they have estimated masses in the narrow range between 1.12 
and 1.18 Mq. Therefore, the question of whether we identify 
a "secondary dip" is legitimate, although strong caution should 
be exercised before claiming it as highly probable. Their mass 



range encompasses four more stars, three of which lie signif- 
icantly closer to the isochrone. If we formulate the intriguing 
hypothesis of an extra-mixing mechanism triggered in this nar- 
row mass range by a low mass companion, we are still left with 
the unexplained case of YBP 1680, which is also likely to have 
a small undetected companion, but does not show any anomaly 
as far as lithium abundance is concerned. This matter deserves 
further attention. 

YBP 890: This star is also a secure member, but, contrary 
to the three discussed above, its distance to the isochrone in the 
CMD is only 0.03 mag. It has a mass of about 1.24 M©, and a 
lithium abundance about 0.9 dex lower than YBP 961, which is 
estimated to be at about the same mass. Its lithium abundance is 
an upper limit, it may therefore be more peculiar than it appears. 
This star has undoubtly undergone a strong amount of extra mix- 
ing. 

YBP 871 and YBP 778 are also significandy under- 
abundant in lithium, and have masses similar to that of YBP 890. 
Their distance to the isochrone in the CMD is about 0.06 mag, 
which is larger than for YBP 890 but not dramatically large. In 
their case, however, we still cannot completely rule out the pos- 
sibility of a combination of errors in the abundance and mass to 
explain their peculiar position in Fig. |2] 

YBP 942: This is a secure member of the cluster, the only 
star with a peculiarly high lithium abundance, which is 0.5 and 
0.8 dex higher than the slightly less massive stars YBP 1320 and 
YBP 1318. Rather than an error in the abundance, an overesti- 
mation of the mass of this star of about 0.05 M© could help us to 
explain its abundance anomaly. However, its distance from the 
isochrone in the CMD is large at about 0.1 mag. 



5.4. Litliium abundance, rotation, and cliromosplieric activity 
in M67 

In stars with a convective envelope, the chromosphere and 
chromospheric activity are generated by non-radiative heat- 
ing m echanisms caused by a magnetic field (see e.g. iHalll 
l2008l for a review). Rotation and differential rotation play a 
key role in this process, because they are at t he base of the 
dynam o effect that sustain magnetic activity (iBohm-Vitensd 
|2007|) . Magnetic braking spins down the star, which there- 
fore loses chromospheric activity as it loses angular momen- 
tum. While until a few years ago there was general consen- 
sus that chromospheric activit y dec ays smo othly during the 
whole stellar life-time ( e.g. iSkumanich 1972; B arrv et anil987l: 
ISoderblom eranil99ll: lDonahuelll998tlLa chaume et al. 19m 
which is still widely accepted dMamaiek & Hillenbrandi2008 ). 
alternative views have also been expressed later on (iPace et al.l 



120091: iLvra & Porto de M ello"2005l: IZhao et al.ll20TTl) . However, 
nobody questions that young stars are more chromos pherically 
active than stars older than 1.5 Gyr (ISoderbloml2010l) . 

This picture suggests that chromospheric activity may be re- 
lated to rotation and age, in a very similar way as lithium abun- 
dance. This is the rationale for a comparison between these three 
parameters in M67. 

The most e xtensive chromo s pheric -activity survey made in 
M67 is that by iGiampapa et al.l (l2006h . who collected data for 
over 60 MS stars. They used the HK index, i.e. a measure of 
the strength of the chromospheric emission in the core of the 
Ca II H and K line, in mA. The HK index is neither trans- 
formed into flux, nor corrected for photospheric contribution. 
Mamaiek & Hillenbrandl ( l2008h transformed the HK indices of 
iGiampapa et aH into R'^^^, i.e. they subtracted the photospheric 
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Fig. 4. On the left: chromospheric activity as a function of mass. On the right: lithium abundances versus chromospheric activity. 
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Fig. 3. Lithium abundance predicted by our models as a function 
of time for different masses. Upper panel: for a MS star with M 
= 1 M0(solid Hne), stars of M = 1.15 M0(dot-dashed) and 1.28 
MoCdotted) currently around the turn-off, and a more massive 
star of 1.33 M0(dashed line) that is already a giant. Lower panel: 
same plot, but after a -0.45 dex shift was applied (as in Fig.|2]i. 



contribution, transformed it into flux, and normalized it for the 
bolometric emission. The data were taken over a time span of 
five years, and measurements for each star were time-averaged 
over large part of this interval, which is of crucial importance 



since the chromospheric activity undergoes a short-term varia- 
tio n analog of the 1 1 -year solar cycle. Thirty-five stars studied 
by iGiampapa et al.l have published lithium abundances that are 
collected here and plotted in Fig. |4] In this figure, we show the 
HK index in mA as a function of mass {Left panel) and HK in- 
dex versus lithium abundances (Right panel). Here upper limits 
are indicated as downward arrows. Circles surrounding the dots 
indicate that the star has a rotation velocity measurement avail- 
able (see below). The left panel has the purpo se of identifying 
the two active outliers. It is similar to Fig. 4 in IGiampapa et all 
in which B-V color is a substitute for mass, but we note that 
among the two active stars in both figures, the only star in com- 
mon is Sanders 1050 /YBP 1342. The star Sanders 747/YBP 681 
is given in Fig. 4 of IGiampapa et al] and not in our Fig. |4] be- 
cause it is a binary, and we excluded it from our compilation. 
On the contrary, we plot San ders 1452/YBP 1067, which is out- 
side the margins of Fig. 4 of IGiampapa et al.[ In the right panel, 
we show that the two active outliers are highly lithium-depleted. 
However, these high levels of lithium depletion are also found 
in many other inactive stars. The Pearson correlation coefficient 
between the HK index and the lithium abundance is weak but 
significant at -0.25. This weak correlation cannot be an effect of 
a correlation between HK index and mass, which, for our sam- 
ple, are completely unrelated. However, the correlation between 
HK index and lithium abundance is entirely due to the two chro- 
mospherically active and highly lithium depleted stars. If we re- 
move them, the correlation between HK and lithium abundance 
is -0.05, which corresponds to a probability of the two quantities 
being unrelated of about 80%. 

iReiners & Giampapal(l2009l) measured the p rojected rotation 
velocities of 15 stars selected from the sample of lGiampapa et al.l 
in order to check whether rotation was at the base of the higher 
chromospheric activity level of some stars. The measurements 
were based on the cross-correlation profile of high-resolution 
spectra. Nine of these stars are present in our compilation of 
lithium abundances, and are indicated in Fig.|4]with a circle sur- 
rounding the dot. With the only exception of YBP 1067, all of 
the projected rotation velocities are consistent with the solar ro- 
tation, as they have an upper limit of about 2 kms"'. On the 
contrary, YBP 1067 has a projected rotation velocity of about 4 
kms"', i.e. it rotates at least twice as fast as the sun. Whether this 
suggests a relationship between rotation and lithium abundance 
is unclear On the one hand, YBP 1067 has a very low lithium 
abundance, on the other hand there are two slow rotators with 
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an even lower lithium content and one more slow rotator with 
an upper limit only slightly above the lithium abundance mea- 
surement for YBP 1067. More data are warranted to shed light 
on this matter. In particular, it would be very important to obtain 
true rotation periods instead of only v sin /. 

6. Conclusions 

Spectroscopic observations of stars in the open cluster M67 pro- 
vide important constraints to test non-standard models of lithium 
depletion. We have rederived the effective temperatures of M67 
stars and corrected the lithium abundances available in the liter- 
ature, in order to have a homogeneous data-set for comparison 
with our models. We have also taken into account NLTE effects, 
providing a homogeneous set of NLTE Li abundances in M67. 
Tables with our T^s, NLTE lithium abundances, and masses are 
presented, so that other groups can test their non-standard stellar 
models. 

M67 stars close to the turn-off, with masses M = L 1 - L28 
Mq, have a peak or a plateau in lithium abundances. Less mas- 
sive stars (M < LI M©) display strong lithium depletion, which 
increases for lower masses owing to a deepening of their con- 
vection zones. Evolved sub-giant and giant stars with M > L28 
Mq show lower Li abundances due to a dilution of their orig- 
inal lithium content. The above pattern is qualitatively well- 
re produced by our models, as well as by the independent models 
of lXiong&Dend (120091) . 

The lithium abundance appears to be a tight function of the 
mass, for stars more massive than the Sun, with a few notable ex- 
ceptions that deserve a closer look. In particular, based on three 
over-depleted stars, we suggest that, at about L15 Mq, the pres- 
ence of a small companion may trigger a large amount of ex- 
tra mixing, or perhaps the initial rotation velocity was higher in 
these stars. More observations are needed before drawing firm 
conclusions about these outliers. 

Our models qualitatively reproduce many observed features 
of the lithium abundance as a function of mass. However, we are 
still far from achieving a close match between observations and 
models with a unique combination of the parameters. 

Two stars in our compilation have chromospheric activity 
levels that are unusually high. They are both highly lithium de- 
pleted but, apart from this circumstance, no strong relationship 
between chromospheric activity and lithium abundance seems 
to be present. Only for one of the two chromospherically active 
stars has a projected rotation velocity been measured, which is 
~ 4 kms"\ the only value that exceeds 2 kms"' among the 
nine measurements in our sample. This suggests that rotation is 
at the base of both lithium depletion and higher chromospheric 
activity in this star, but that other extra mixing processes must 
also be efficient in some slow rotators. 
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Table 1. M67 data. See Sects. 2 and 3 for explanations. 
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5752 
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99 
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0.1 
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14.358 
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97 
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14.947 
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5495 
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5970 
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99 


33.206 
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1.7 
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0.1 
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14.559 
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34.67 


3 


1.0 


5929 
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5727 


0.94 
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0.1 


YBP 1075 


13.712 


0.674 


99 


32.974 


3 


0.7 


5871 


0.94 


6123 


0.69 


5749 


0.1 


YBP 1088 


14.492 


0.659 


93 


33.199 


3 


1.0 


5890 


0.91 
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0.93 
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up.lim. 
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V 
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T 
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(mag) 


(mag) 


% 
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(dex) 


(K) 


(dex) 


(K) 
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(K) 


dex 


YBP 1090 


13.8 


0.65 


100 


32.333 


3 


2.3 


6086 


2.31 


6095 


2.11 


5824 


0.1 


YBP1129 


14.171 


0.624 


99 


34.528 


3 


2.1 


5959 


2.09 


5942 


2.06 


5907 


0.1 


YBP 1137 


14.873 


0.698 


97 


33.851 


3 


1.0 


5741 


0.87 


5546 


0.95 


5675 


up.lim. 


YBP 1197 


13.315 


0.606 


100 


34.526 


3 


2.2 


6207 


2.14 


6123 


2.03 


5967 


0.1 


YBP 1247 


14.144 


0.609 


99 


32.991 


3 


2.1 


5994 


2.06 


5942 


2.07 


5957 


0.1 


YBP 1334 


14.403 


0.68 


99 


32.923 


3 


2.0 


5957 


1.89 


5807 


1.82 


5730 


0.1 


YBP 1387 


14.098 


0.626 


99 


33.828 


3 


2.3 


6090 


2.22 


5970 


2.16 


5901 


0.1 


YBP 1458 


14.977 


0.739 


98 


33.351 


3 


1.0 


5640 


0.87 


5482 


0.94 


5555 


up.lim. 


YBP 1496 


13.879 


0.607 


100 


34.607 


3 


2.6 


6173 


2.52 


6053 


2.45 


5963 


0.1 


YBP 1504 


14.171 


0.625 


99 


33.733 


3 


2.2 


5934 


2.21 


5942 


2.18 


5904 


0.1 


YBP 1514 


14.777 


0.721 


98 


33.878 


3 


1.6 


5613 


1.59 


5597 


1.6 


5607 


0.1 


YBP 1587 


14.163 


0.641 


99 


32.336 


3 


2.0 


5975 


1.98 


5942 


1.91 


5852 


0.1 


YBP 1622 


14.156 


0.632 


99 


33.407 


3 


2.2 


6043 


2.13 


5942 


2.08 


5881 


0.1 


YBP 1716 


13.299 


0.619 


96 


33.948 


3 


2.3 


6030 


2.36 


6123 


2.22 


5924 


0.1 


YBP 1722 


14.13 


0.601 


96 


34.007 


3 


2.2 


6007 


2.16 


5956 


2.18 


5983 


0.1 


YBP 1735 


14.332 


0.661 


97 


33.104 


3 


0.9 


5959 


0.83 


5847 


0.81 


5789 


0.1 


YBP 1758 


13.207 


0.653 


98 


34.398 


3 


2.4 


6221 


2.29 


6067 


2.11 


5814 


0.1 


YBP 1768 


14.404 


0.656 


99 


34.222 


3 


2.1 


5844 


2.07 


5807 


2.07 


5805 


0.1 


YBP 1788 


14.441 


0.663 


84 


33.763 


3 


1.8 


5886 


1.72 


5780 


1.72 


5783 


0.1 


YBP 1852 


13.962 


0.613 


99 


32.324 


3 


2.1 


6009 


2.11 


6025 


2.05 


5943 


0.1 


YBP 1903 


14.733 


0.689 


99 


32.409 


3 


1.0 


5609 


1.01 


5623 


1.07 


5703 


up.lim. 


YBP 1948 


14.015 


0.612 


97 


33.136 


3 


2.4 


6164 


2.29 


5997 


2.25 


5947 


0.1 


YBP 1955 


14.212 


0.63 


98 


32.588 


3 


2.2 


5961 


2.17 


5915 


2.14 


5888 


0.1 


YBP 901 


13.422 


0.554 


100 


- 


4 


2.63 


6191 


2.57 


6137 


2.58 


6145 


0.07 


YBP 963 


12.76 


0.565 


100 


- 


4 


2.1 


6210 


2.09 


6194 


2.03 


6106 


0.08 


Sand. 998 


13.06 


0.558 


- 


- 


4 


2.64 


6223 


2.5 


6053 


2.56 


6131 


0.06 


YBP 871 


13.152 


0.582 


100 


- 


4 


1.86 


6156 


1.79 


6053 


1.79 


6048 


0.16 


Sand. 1064 


14.04 


0.661 


- 


32.7 


5 


2.14 


5845 


2.28 


5997 


2.12 


5789 


0.1 


YBP 713 


14.172 


0.714 


99 


35.443 


5 


2.07 


5800 


2.2 


5929 


1.97 


5628 


0.1 


YBP 1105 


13.946 


0.59 


99 


32.8 


5 


2.55 


6094 


2.5 


6039 


2.49 


6020 


0.1 


YBP 1397 


14.009 


0.622 


97 


36.808 


5 


2.28 


5972 


2.32 


6011 


2.25 


5914 


0.1 


YBP 1342 


14.285 


0.65 


99 


33.684 


5 


1.15 


5810 


1.24 


5874 


1.2 


5824 


up.lim. 


Sand. 1057 


14.3 


0.668 


- 


34.5 


5 


1.75 


5818 


1.82 


5861 


1.75 


5767 


0.1 


YBP 1486 


13.864 


0.574 


94 


- 


6 


2.63 


6178 


2.54 


6067 


2.54 


6075 


0.1 


YBP 890 


13.179 


0.59 


100 


34.212 


6 


1.48 


6153 


1.43 


6067 


1.4 


6020 


up.lim. 


YBP 961 


13.19 


0.576 


100 


34.59 


6 


2.37 


6151 


2.31 


6067 


2.31 


6068 


0.1 


YBP 1680 


13.646 


0.645 


100 


35.428 


6 


2.45 


5934 


2.6 


6137 


2.38 


5840 


0.1 


Sand. 1055 


13.8 


0.593 


- 


- 


6 


2.43 


6116 


2.41 


6095 


2.35 


6010 


0.1 


YBP 1456 


12.841 


0.943 


100 


- 


7 


0.5 


5000 


0.75 


5105 


0.67 


5008 


up.lim. 


YBP 877 


12.732 


0.583 


99 


- 


7 


1.9 


6205 


1.89 


6194 


1.79 


6044 


up.lim. 


YBP 891 


11.401 


1.089 


98 


- 


7 


0.0 


4740 


0.2 


4731 


0.18 


4710 


up.lim. 


YBP 1017 


12.361 


0.962 


96 


33.639 


7 


0.5 


4905 


0.7 


4931 


0.74 


4962 


up.lim. 


YBP 1024 


9.586 


1.443 


98 


31.655 


7 


-1.0 


4010 


-1.01 


4295 


-0.83 


4176 


up.lim. 


YBP 1337 


12.878 


0.972 


100 


- 


7 


0.9 


4850 


1.35 


5105 


1.18 


4940 


0.1 
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Table 2. Hyades data 
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vB 


V 


B-V 


Teff 


logg 


Mass 


EW 


AEW 


A(Li) 


AA 


EWref 


phot ref 


2 


7.78 


0.617 


5916 


4.36 


1.140 


89.0 


2 


2.73 


0.03 


T 


Bal. 


3 


8.88 


0.848 


5099 


4.48 


0.847 


<25.0 




<1.37 




D 


Bal. 


4 


5.97 


0.341 


6897 


3.99 


1.524 


<3.0 




<1.79 




T 


Bal. 


6 


5.95 


0.348 


6861 


3.96 


1.513 


54.8 


3 


3.10 


0.04 


BB 


Joner 


8 


6.37 


0.419 


6544 


4.02 


1.411 


<4.0 




< 1.69 




BB 


Bal. 


10 


7.85 


0.589 


6000 


4.42 


1.182 


82.0 


2 


2.74 


0.03 


T 


Bal. 


13 


6.62 


0.420 


6542 


4.12 


1.410 


<5.0 




<1.79 




BT 


Bal. 


14 


5.73 


0.355 


6825 


3.87 


1.502 


65.0 


3 


3.18 


0.04 


BT 


Bal. 


15 


8.09 


0.658 


5703 


4.39 


1.037 


57.0 


2 


2.32 


0.03 


T 


Bal. 


17 


8.46 


0.696 


5539 


4.47 


0.971 


42.0 


2 


2.03 


0.05 


T 


Bal. 


18 


8.06 


0.638 


5779 


4.41 


1.072 


74.0 


2 


2.52 


0.04 


T 


Bal. 


19 


7.14 


0.512 


6231 


4.22 


1.290 


82.4 


3 


2.91 


0.03 


BB 


Bal. 


20 


6.32 


0.399 


6624 


4.03 


1.438 


87.6 


3 


3.21 


0.03 


BB 


Bal. 


21 


9.15 


0.816 


5207 


4.63 


0.872 


<3.0 




<0.58 




Ra 


Bal. 


26 


8.63 


0.743 


5437 


4.51 


0.936 


13.0 


2 


1.39 


0.07 


T 


Bal. 


27 


8.46 


0.715 


5506 


4.46 


0.960 


23.0 


2 


1.72 


0.04 


T 


Bal. 


31 


7.47 


0.566 


6067 


4.29 


1.215 


95.0 


2 


2.88 


0.02 


T 


Bal. 


36 


6.80 


0.441 


6463 


4.17 


1.382 


5.9 


3 


1.81 


0.18 


BT 


Bal. 


37 


6.61 


0.405 


6600 


4.14 


1.430 


10.3 


3 


2.14 


0.12 


BT 


Bal. 


38 


5.72 


0.320 


7019 


3.93 


1.560 


19.6 


3 


2.70 


0.14 


BT 


Bal. 


42 


8.86 


0.759 


5414 


4.59 


0.929 


11.0 


2 


1.30 


0.08 


T 


Bal. 


44 


7.19 


0.450 


6433 


4.31 


1.371 


20.4 


3 


2.34 


0.07 


BT 


Bal. 


46 


9.11 


0.867 


5028 


4.55 


0.832 


<4.0 




<0.55 




T 


Bal. 


48 


7.14 


0.521 


6203 


4.21 


1.278 


91.0 


2 


2.95 


0.02 


T 


Bal. 


51 


6.97 


0.443 


6457 


4.23 


1.380 


6.5 


3 


1.84 


0.17 


BT 


Bal. 


49 


8.24 


0.585 


6012 


4.58 


1.188 


57.0 


2 


2.55 


0.03 


T 


Bal. 


59 


7.49 


0.543 


6136 


4.33 


1.248 


82.0 


2 


2.84 


0.03 


T 


Bal. 


61 


7.38 


0.514 


6224 


4.32 


1.287 


110.8 


3 


3.08 


0.03 


BB 


Bal. 


64 


8.12 


0.657 


5707 


4.40 


1.039 


59.0 


2 


2.34 


0.03 


T 


Bal. 


65 


7.42 


0.535 


6160 


4.31 


1.259 


106.0 


2 


3.00 


0.03 


T 


Bal. 


66 


7.51 


0.555 


6100 


4.32 


1.231 


73.0 


2 


2.75 


0.03 


T 


Bal. 


73 


7.84 


0.609 


5940 


4.39 


1.152 


85.0 


2 


2.72 


0.03 


T 


Bal. 


77 


7.05 


0.500 


6263 


4.20 


1.304 


19.0 


6 


2.20 


0.13 


Re 


Bal. 


76 


9.20 


0.759 


5392 


4.72 


0.922 


16.0 


2 


1.45 


0.06 


T 


Bal. 


78 


6.92 


0.453 


6422 


4.20 


1.367 


32.6 


3 


2.56 


0.05 


BB 


Bal. 


79 


8.93 


0.827 


5171 


4.53 


0.863 


<3.0 




<0.56 




T 


Joner 


81 


7.10 


0.470 


6365 


4.25 


1.345 


15.6 


3 


2.18 


0.08 


BT 


Bal. 


85 


6.51 


0.426 


6519 


4.07 


1.402 


<6.0 




<1.85 




BB 


Bal. 


86 


7.05 


0.463 


6388 


4.24 


1.354 


20.7 


3 


2.32 


0.07 


BT 


Bal. 


87 


8.58 


0.743 


5437 


4.49 


0.936 


12.0 


2 


1.36 


0.07 


T 


Bal. 


88 


7.75 


0.554 


6103 


4.42 


1.232 


89.0 


10 


2.86 


0.06 


S 


Joner 


90 


6.40 


0.413 


6568 


4.05 


1.419 


<3.0 




< 1.58 




BB 


Bal. 


92 


8.66 


0.741 


5443 


4.52 


0.938 


15.0 


2 


1.46 


0.06 


T 


Bal. 


93 


9.40 


0.883 


4968 


4.64 


0.820 


<3.0 




<0.29 




T 


Bal. 


94 


6.62 


0.431 


6499 


4.11 


1.395 


<2.2 




<1.40 




BT 


Bal. 


97 


7.93 


0.634 


5793 


4.36 


1.079 


84.0 


2 


2.60 


0.05 


T 


Bal. 


99 


9.38 


0.851 


5090 


4.68 


0.845 


4.0 


2 


0.59 


0.13 


T 


Bal. 


101 


6.65 


0.433 


6493 


4.12 


1.393 


<3.0 




<1.53 




BB 


Bal. 


105 


7.53 


0.575 


6042 


4.31 


1.203 


87.0 


2 


2.81 


0.02 


T 


Bal. 


109 


9.40 


0.817 


5203 


4.73 


0.871 


5.0 


1 


0.74 


0.09 


Ra 


Bal. 


116 


9.01 


0.821 


5191 


4.57 


0.868 


<5.0 




<0.73 




T 


Bal. 


118 


7.74 


0.580 


6026 


4.39 


1.195 


77.0 


2 


2.72 


0.03 


T 


Bal. 


127 


8.92 


0.710 


5520 


4.65 


0.964 


18.0 


2 


1.62 


0.05 


T 


SIMBAD 


153 


8.91 


0.859 


5062 


4.48 


0.839 


8.0 


2 


0.82 


0.10 


T 


Bal. 


121 


7.29 


0.504 


6256 


4.29 


1.301 


114.9 


3 


3.12 


0.02 


BT 


Bal. 


124 


6.25 


0.501 


6261 


3.88 


1.303 


9.0 


3 


1.86 


0.13 


BT 


Joner 


128 


6.75 


0.450 


6433 


4.14 


1.371 


14.0 


3 


2.17 


0.09 


BT 


Bal. 


180 


9.10 


0.853 


5081 


4.57 


0.843 


5.0 


2 


0.65 


0.13 


T 


Bal. 


187 


8.60 


0.776 


5352 


4.46 


0.912 


17.0 


1 


1.44 


0.04 


Ra 


Joner 
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